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SUMMARY 


In this report calculations of the ratios of the velocity increments 
at test bodies to those at the tunnel walls caused by the solid blocking 
of these bodies within the walls of a closed rectangular wind tunnel are 
presented. The bodies treated include two-dimensional airfoils 5 small 
bodies of revolution] straight, untapered, finite-span wings of varying 
span] and swept, untapered, finite-span wings of varying span. It is 
shown that, after wake blocking effects have been removed, the present 
method furnishes semiempirical blocking corrections for most wind-tunnel 
models and their components. Res-ults are presented for all the cases 
mentioned. The test-section proportions of the Southern California 
Cooperative Wind Tunnel at the California Institute of Technology ^viz., 
ratio of height to width equal to l/\[ 2 ) are used in calculations. 

INTRODUCTION 


The velocity-ratio method of obtaining blocking corrections in 
high-speed, subsonic wind tunnels was first solved by Gdthert for the 
cases of a body of revolution and of a finite-span wing of span-to- 
diameter ratios of 0.25 and 0.5 in a closed circular tunnel (reference 1 ) 
This work was later extended in an unpublished report to the cases of a 
wing having a span-to-diameter ratio of 0.?5 and of a wing spanning a 
closed circular wind tunnel. 

In the present report, the methods of reference 2 are used to extend 
the previous results for straight wings of varying span in a closed rec- 
tangular tunnel whose height-to-width ratio is For swept wings, a 

slightly different approach involving the use of line doublets has been 
utilized. 

This work was done by the California Institute of Technology and 
has been made available to the NACA for publication because of its 
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general interest. The author wishes to express his thanks for the 
assistance rendered by Dr. Clark B. MiUikan, Mr. Josiah E. Smith, 
and Mr. Richard W. Bell in the preparation of this report, and also to 
thank Misses Dorothy Lodter and Donna Deeney for their performance of 
the extensive computations. 


SYMBOLS 


a 

A 

c 

d 

g 

H 

I 

M 

Mu 


m 

n 

q 

q' 

r 


R 

s 

s' = s/W 


distance of point source and sink from doublet origin 
center, or axis of test section 
wing chord 

maximum diameter of a body of revolution 
x-coordinate of a source segment 
tunnel height 

length of a body of revolution 
Mach number, corrected for blocking 

tunnel calibration Mach nximber, uncorrected for blocking 
(may include strut calibration) 

integer, indicating image nvimber in y-direction 

integer, indicating image number in z-direction 

line-source strength, square feet per second 

point-source strength, cubic feet per second 

distance from source or doublet element to point at which 
velocity increment shall be obtainedj j -2 s + y^ + 

remainder 

semispan of model wing (measured in y-direction) 


t 


maximum wing thickness 
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u total axial velocity increment in test section due to all 

images except primary one (model) 

u' total axial velocity increment in test section due to 1 

images, including primary one (model) 

U]^ axial velocity increment in test section due to a single 

image 

U tunnel axial velocity 

v^ velocity in x-direction due to a point doublet 

W tunnel width 

w radial coordinate; = y^ + 

X axial coordinate 

X]_ = X 

y lateral coordinate 

y^ = y - mW 

z ■ vertical coordinate 

z^ = z - nH 

e = sin“^ ^ 
r 

M. line-doublet strength, cubic feet per second 

p-' point-doublet strength, feet^ per second 

A angle of sweep of a given wing at any Mach number My^ 

Aq equivalent angle of sweep at My = 0 

■'I'’ Stokes stream function, cubic feet per second 
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Subscripts : 

A, B, C, and particular points in test section or on test-section 

so forth walls 

For further explanation of the symbols see figures 1 to U* 


DETERMINATION OF INTERFERENCE VELOCITIES 
Two-Dimensional Wing 


A two-dimensional wing may be represented by a chordwise distri- 
bution of infinite line sources and sinks. The axial velocity increment 
produced by any single infinite-line-source image (fig. 1) is 


^ sin e = 

1 2nr 2itr2 


( 1 ) 


at a point A in the center of the basic tunnel, r*^ = n^H^ + g^. Thus, 


u 


qg/ 1 


qg / 1 


^A 2Tt(^2jj2 + g2y 2nH^| 2 ^ 


( 2 ) 


Omitting the central source and summing for the remainder gives 


_ qg 

p / p 

n=l 1 „2 ^ gf 

h2- 


(3) 


Similarly, the axial velocity at the wall of the tunnel, point B, 
including both the central source and all the images, is 


u 


_ qg 


I 


nH n=>l 



(U) 
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For normal chord 
and hence 
with the result 


sizes, g is small compared with the tunnel dimensions 
may be neglected in comparison with n^ or ^n - , 



( 5 ) 


Thom has shown in reference 2 that qg for a single line source may be 
replaced by ^qg, the distribution of sources and sinks representing 
the wing section. Therefore, since the ^qg terms would also cancel 
each other, the result obtained for a single line source in equation (5) 
is identical to that for a complete wing. 


Body of Revolution 


A body of revolution may be represented by a distribution of point 
sources and sinks along the tunnel center line. The axial velocity 
increment due to a single image point source (fig. 2) is 


q' 


sin 


e = 


ij.nr3 


(6) 


? ? ? ? ? ? 

at point A in the center of the basic tunnel, = g“^ + n^H + 

Then, substituting into equation (6), 


U-, 


= q'g 


/ii . n2 4 - -2} 

\W^ / 


\3/2 


(7) 


Again omitting the central source and summing for the remainder with the 

2 

assumption that g2/w2 is negligible compared with n2 2^ 


and m2 gives 
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^A = 


q'g 


iinW^ 


00 00 

uZZ/ 


,3/2 


+ 2 


n=l m=l U2 SI H- m2 


00 

z 


1 \- 


n=l n2 ll 1 ' 

(m=0)y (n=0) 


3/2 

s'! * “ fei 


00 . V 

Z ^ 


3/2 


\ W“= 


( 8 ) 


Similarly, it can be shown that the wall velocities u'g and u'q are 


u' 




00 00 
fZZ 

n=l m=l 

1 

3/2 ~ 

" 2 2_ 
n=l 

1 

3/2 

L . if . ,.2 

(n - if £ 

• f 


\ 2 ; w2 J 

(m=0) 




( 9 ) 


and 


u'p = 


_ q' 


linW^ 


00 00 

ss 


n 


R / 1\ 


3/2 ^ 

. 22_ 

m=l 


(n=0) 


-l3/2 


(m-l 


> 


(10) 


The latter two equations include the effect of the primary source, which 
must be omitted in the calculation of u^. 


Straight, Untapered, Finite-Span Wing 

A finite-span wing may be represented by a distribution of finite- 
length soiirces and sinks. The axial velocity increment produced by a 
single source element is (fig. 3), for this case. 


du^ 


q dy^ 
Unr^ 


sin 9 


qg dy^ 


(11) 


For the point A in the center of the basic tunnel^ the general expres- 
sion for the square of the distance from the source element is 

r 2 ^ 

the total contribution at A of one images namely^ 


n2n2 


^mW + Integration across the image span gives 
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ui =p 
-^A kn , 




f-S 


J^g2 + + ^mW + 

Performance of the integration leads to the result 


3/2 


( 12 ) 


Un 


qg 

mW + s 

lin(n^H^ +g^) 

[g2 + + (mW + s)^]]^'^^ 


mW - s 


I 


J 

(13) 


which is the same as equation I 6 , reference 2. Making equation (13) 

g2 

nondimens ional and again neglecting the terms give 



( 111 ) 


As before for the case of the body of revolution, the total velocity 
increment is 


^A 


00 00 


n=l m=l m=l n=l 

(n-O) (m-0)J 




(13) 


However^ it will be noted that the single summation for which n = 0 
leads only to the indeterminate quantity This may be evaluated by 

the application of Hospital *s Rule^ finally giving for a single 
source line 
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Similarly, the total axial velocity increments at points B and C on 
the walls of the tunnel are 
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and 



( 18 ) 

Furthermore, at ar^r point y along the span of the wing, the result is 



As before, the primary source line which corresponds to the wing is not 
included in calculating the velocity increments at any point in the wing. 
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Swept, Untapered, Finite-Span Wing 

In the preceding treatment for unswept bodies the analyses were 
carried out for single line and point sources as a simplification for 
bodies which could be represented by a combination of line and point 
sources and sinks. Thus, for a given value of g, all the image 
sources are in the same plane and a constant distance from the plane 
of A, B, and G, namely, g. Obviously, this condition no longer holds 
for swept bodies in which g would be a function both of span and 
angle of sweep, and hence the simplified treatment is no longer appli- 
cable, It is then necessary to proceed to a more general representation 
of the flow field involving the use of both sources and sinks. One of 
the simplest of such combinations is the doublet and that is what is 
used. 


The Stokes stream function for a point doublet is given by 
(reference 3 ) 


t =-M-' 



( 20 ) 


where . The point-doublet strength p* remains finite as a. 

kn 

the distance of the point source and sink from the doublet origin^ 
approaches zero and q* goes to infinity* The velocity in the x-direction 
caused by the doublet element^ whose axis is parallel to the x-axis 
(fig. ij.), is given by 



( 21 ) 


Then, for a swept-wing element, with doublet-element axes parallel to 
the x-axis. 


du^ = 



( 22 ) 
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The square of the distance from the doublet element to A is 

+ ^mW + Considering Aq positive for sweepback, 

then for the right wing with respect to the model, — = tan Aq, and 

-X2_ ^1 

for the left wing, = tan Aq. Having eliminated x^_, the expres- 

^1 

sion for the contribution of a single-image doublet is 




n^H^ + ^mW + 


y^^tan^AQ + n^H^ + ^mW + yq^ ^ 


i5/2 


> dy. 


(23) 


It should be noted that a change of sign between the right and left 
halves of a wing makes no difference mathematically. Therefore, the 
swept wing acts in exactly the same manner as a wing yawed at an 
angle Performing the integration and letting 


= n2 SI . m2 

W2 


b^ ' = 2m 

p 

c = tan Aq + 1 


dA' = Sa^'c - (b^')' 


s ' = 


W 


%+' = ^a' ^ ^a'®' ^ 
V = ^a' - ^a'2' + 
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The total velocity increment at A may be obtained ty the same summation 
as indicated for the straight wing in equation ( 13 ). However, for m 
and n equal to zero, the constants given above become; 




and, for n = 0, 
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The final result is then 
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In following a similar procecure to obtain u'g and u'q, it is found 
that the constants are the same as for except that the subscripts 

are changed to B and C, respectively, and n — > ^ ^ ® 

m — > m - ^ at C. The summations are the same as for the straight 
wingj namely, 
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and 



Support Struts 

The method previously used to calculate the velocity ratios for 
aerodynamic bodies may also be applied to support struts. However, in 
order to avoid infinite velocity increments at the junctions of struts 
and model, it is necessary to consider the support system in the tunnel 
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as an integral part of the model. With this viewpoint, velocity ratios 
could also be calculated for struts in the same manner as for wings, in 
which just the images are summed in order to calculate Uj^. The span 

of the struts would be constant at 2s/W =1.0 in a reoriented tunnel 
for which H/W = \[2. 

Using points A and C, where C would now be at the top of the reori- 
ented tunnel, the velocity ratio %/u'q for half wings or struts would 

be identical to those for complete wings or double struts, both com- 
pletely spanning the tunnel. There is one simple case for which the 
result is immediately known, namely, a single, centrally mo\mted, 

Tinswept strut for which the two-dimensional result of l/3 applies (see 
section Two-Dimensional Wing). In general, if other than a single-strut 
support system were used, the velocity ratio would be a function of the 
strut spacing used as well as of the angle of sweep. The problem would 
be further complicated by the presence of a rear strut, which is fre- 
quently the case. 

Because of the additional complexity involved and the expected 
difficulty in separating the total wall velocity increment into the 
separate effects due to the solid blocking of the support system, model 
wing, and model fuselage and to the wake blocking of each component, all 
of which may have different velocity ratios, no general solution of the 
■support-strut case has been presented. It would be sin?)ler and probably 
more acciirate to perform a complete calibration of the area in which a 
model would normally be mounted, with the struts installed. During this 
calibration, the wall pressures at B and C could also be obtained, thus 
giving base values which include both the solid and wake blocking and 
interference effects of the model support system. 


NUMERICAL CALCULATIONS AND RESULTS 


The methods used in slimming the doubly infinite series are 
eiqDlained in detail in appendix A. Briefly, calculations are made for 
each image up to a finite number n = m = n^^. The remainders are 

obtained by direct integration from n^^ to infinity after making certain 

simplifying assumptions. Unfortunately, the series convergence is not 
very rapid and it is necessary to take n^^ as high as seven in most 

cases. Fortunately, however, there is a negligible difference between 
the remainder terms for the swept and straight wings, since the effect 
of sweep rapidly diminishes as the distance between images and tunnel 
increases. Hence the very difficult problem of attempting to integrate 
the complicated remainders for the swept wings is avoided. In all 
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cases, the dimensional ratio of the Southern California Cooperative 
Wind Tunnel, namely, H/W = l/\/^, has been used. 

Figures 3 and 6 show the variation with 2s/W of the velocity 
ratios u-a/u.'b and Ua/'u'q, respectively, for sweep angles of 0°, 30°, 

U5°, and 60°. It should be noted that each set of curves begins from a 
common value at 2s/W = 0, which corresponds to the body- of -re volution 
case, and that the effect of swee^ is to increase the ratios over the 
values for Aq = 0. The effect of the side walls upon ua/u'b and 

ua/u'c for straight wings as 2s/W approaches unity is also of 
some interest. In. the former case, as the wing approaches the side 
walls their effect upon continues until the wing just touches 

the walls, at which point the ratio discontinuously drops from the 
three-dimensional result which includes the effect of the side walls to 
the purely two-dimensional value of 1/3 in which the side-wall effect 
has been eliminated. As the wing tip nears the side wall at point C 
(2s/W — > 1 ) it also causes u'q to approach infinity while u^ remains 

finite, so that u^/u'q approaches zero. 


It should be noted that, although the velocity ratios for swept 
wings are greater than those for straight wings, the absolute magnitudes 
of Vlj^ and u'g decrease with increasing sweep for a given span. The 

value of u'q is also less for Aq / 0 than for Aq = 0, but behaves 


irregularly beyond 2s/W = 0.7. These three blocking velocity incre- 


ments, made nondimensional by multiplication with 

in figures 7, 8, and 9 as a f-unction of wing span, 
given for these wings are based upon line doublets. 


1 ^ 


are plotted 


The numerical values 
The corresponding 


values of u'g 


and 



for a body of revolution. 


which is represented by a point doublet, are 16.2, 36.8, and 29. i;, 
respectively. From figure 9 it can be seen that u'q becomes negative 


for small angles of sweep as the wing tip approaches the wall 
(2s/W-^ 1.0). 


In reference 2, Thom has used an arbitrary line or point source of 
strength q displaced a distance g from the origin, giving the 
multiple qg which appears in his equations. Thom considered two 
elementary shapes in two and three dimensions, namely, a Rankine oval 
and ovoid, and the shapes corresponding to a line and point source of 
strength +q at a distance g from the origin, a line and point sink of 


strength 



at the origin, and a line and point sink of strength 
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at a distance -g from the origin. For the Rankine oval and ovoid^ 
which are the bluntest bodies obtainable by this method, ~3~qg is 2qg, 
which must be a maximum. For the latter cases, ^qg = ^ qg- 

/ Q 

calculated g/c and as a function of thickness ratio t/c for 

g/l and 

Ud^ 


TO 

both two-dimensional shapes and 


Thom has 


as a function of d/l 


for both body-of-revolution shapes. This procedure may be repeated for 
any arbitrary soixrce-and-sink distribution to facilitate a more correct 
calculation of u'g/U, and 'u'q/U by Thom's method for a given 

shape of wing or fuselage. 


Although the numerical values of figures 7, 8, and 9 are based 
upon line doublets, they have been divided by ti. and hence are equiva- 
lent to the results for a line source divided by qg, since p = 2qg. 
Therefore Thom's method, wherein ^qg for a specific shape replaces 
qg or p when using these figures, is still directly applicable. 

Using the values of and g/c given in figure U of reference 2 

for a Rankine oval, one can easily calculate the incompressible-flow 
values of u^^/U, u'g/U, and u'q/U. For example, for a straight wing 

of this section having unit chord, thickness ratio t/c = 0.1, and 
span 2s/W = 0.8 in a tunnel of width W = 12 feet, Ua /U = 0.00075. 


APPLICATION OF METHOD 


In the preceding sections, ratios of the velocities at the model 
to those at reference points on the wall have been calculated for 
certain representative bodies. During a given test, then, the velocities 
at the walls must be measured in the model plane ABC, and so forth, and 
compared with the reference values obtained without the model installed. 
The difference of the two is not immediately u'g or u'q, however. 

Included in the measurement are the effects of the struts, if the model 
is held in this manner and if the base wall pressvires do not include 
the strut effect, and also an amount due to wake blocking. It will be 
assvimed that the tunnel has been calibrated with the struts installed. 

The wake blocking presents a more complicated problem, but its magnitude 
is fortunately small compared with the solid blocking for most models. 

The methods of separating the wake and solid blocking effects and 
obtaining the final corrected dynamic pressure and Mach number M from 
the measured wall pressures have been outlined in reference U. The 
procedure given therein for the separation of wake and solid blocking 
may be somewhat modified for straight wings and fuselages by the use of 
the wall pressures rather than use of the theoretical wake blocking as 
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given by Thom (reference 2), The latter has represented the wake of a 
body of revolution or a wing by a point or line source, respectively, 
at the body with a corresponding point or line sink downstream at 
infinity. If the sources are assumed to be in the model plane, their 
axial contribution in that plane is zero. This constitutes a good 
assumption for a body of revolution and a straight wing. Thus for these 
cases the only wake effects at the model plane are. due to the sinks at 
infinity, which create a \iniform disturbance in the flow field. There- 
fore, the velocity ratios due to wake blocking are unity for straight 
wings and three-dimensional bodies. The absolute magnitude of the wake 
blocking effect may be obtained by the fact that the solid blocking 
effects at the walls disappear downstream of the model and the remaining 
velocity increment approaches the total disturbance due to the point or 
line sources and sinks representing the wake. The wall wake increment 
at the model plane is then just one-half of the total far downstream 
because of the zero contribution of the sources at the model plane. 
Therefore, the wall wake increment may be determined and subtracted from 
the total velocity increments at B and C, enabling the solid blocking 
corrections to be evaluated. 


The experimental procedure outlined above is not strictly appli- 
cable to the case of swept wings because they extend downstream from 
the reference plane ABC, which only passes through the apex of the 
wings. Therefore, the velocity ratios due to the wake will no longer 
be unity because of the finite contribution made by the wing source 
elements to the wall velocities at B and C. However, for wings of small 
sweep and of small span compared with the tiuonel width combined with 
fuselages, for which the previous argument still holds, it may be 
expected that the velocity ratios due to wake blocking in the apex 
plane ABC are sufficiently close to unity to permit use of this procedure 
with sufficient accTiracy. 


After the strut and wake effects upon the wall velocity have. been 
eliminated, the velocity ratios given in this report may be applied. 
However, the factors for wing and fuselage (body of revolution) are in 
general dissimilar, and it is necessary to estimate the relative magni- 
tude of the correction due to each. This may be done from the model 
solid blocking results of reference 2 at the model for fuselages and 
straight wings or by comparing the results of figure 7 with the body- 


of -revolution result of 



17.0, knowing ^qg for both wing 


and fuselage 5 the estimated relative magnitudes at the wall may be 
obtained by application of the respective velocity ratios. As mentioned 
previously, sweep decreases the velocity increments at all points in the 
flow field. Figure 10 shows the variation of the ratio of u^ with 


sweep to Vi^ for Aq = 0 as a function of Aq for various spans. 
This figure will facilitate estimation of the relative magnitudes of 
the fuselage and wing corrections when Aq ^ o°. 
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Now that the meastired values of u'/U due to solid blocking at 
points B and C on the walls have been split up into the separate contri- 
butions of the wing and fuselage, u/U at point A in the model for 
each is known. However, the increase in velocity is generally of lesser 
interest than the corresponding increase of Mach number and dynamic 
pressure. For small corrections, the equations given by Thom (refer- 
ence 2) and also in reference k, revised to the present notation, are 



The ratio q/q^ is now the ratio of the c^mamic pressure corrected for 
solid blocking to the uncorrected, or calibration dynamic pressure. 

Using the above corrections and a calculated value of the wake 
blocking for a given model, it is possible to operate at a desired 
corrected Mach number at the model by setting a predetermined wall Mach 
number at point B or C. This procedure is particularly desirable in 
tunnels where rapid determination of the test results is required, such 
as at the CWT. 

The velocity-ratio method of obtaining blocking corrections has 
certain advantages over previously used methods. First of all, being 
semiempirical in nature, it promises greater accuracy, particularly 
since the velocity ratios are usually less than unity. Hence the cor- 
rected Mach numbers and dynamic pressures are less than the measured 
ones, whereas by other methods they are higher, being the result of 
increases of lower observed values by theoretical correction factors. 
Procedures which raise observed Mach numbers to higher corrected values 
by means of theoretical formulas are always inaccurate and dangerous at 
subsonic speeds closely approaching M = 1. However, correcting down 
from a high observed Mach number to a lower corrected value as is done 
with the wall-pressure procedure is much safer and more accurate. 
Secondly, as Gothert has pointed out in reference 1, the Prandtl-Glauert 
rule no longer applies at the body in the supercritical speed range but 
is again applicable some distance from the model. In fact, Gothert 
states that the velocity-ratio method may be applied with good accuracy 
up to the speed at which the shock waves approach the tunnel wall, 
namely, the choking speed. Most blocking theories use the Prandtl- 
Glauert rule and, of course, break down as the critical Mach number of 
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the model is exceeded. Furthermore, although the present analysis has 
been carried out at zero lift, Gbthert has pointed out in reference 1 
that the procedure is applicable to a first approximation to cases with 
lift by averaging the pressures at the two walls opposite the lifting 
surfaces. This would normally be at points B and E. 

The derivations of the preceding sections have been carried out 
using line and point sources and doublets. These unit singularities 
may be replaced by combinations of sources and sinks properly dis- 
tributed to represent the desired shape. The factors qg and p 
determine the body shape for straight wings and bodies of revolution 
and for swept wings, respectively. Replacement of these factors by 
^qg and ^p, which are chosen to simulate the desired shape, leads 
to the actual velocity increments in the incompressible-flow field for 
this particular body. However, it is seen that these shape factors are 
common multiples in all equations for each case 5 therefore the velocity 
ratios are independent of shape. 

It is also shown in reference 1 that a given size of compressible- 
flow field corresponds to a fictitious incompressible one in which 
the axial (x) dimension is unchanged and the lateral and vertical 

(y and z ) dimensions are shrunk by an amount y^ = \jl - and 

Zi = Zq \|l - where the subscripts c and i refer to com- 

pressible and incompressible, respectively. Therefore, for straight 
wings and bodies of revolution, ^qg and ^p are functions of Mach 
number, but the velocity ratios in the plane of the origin (x = O) 
remain constant. This also applies to swept wings as far as airfoil 
shape is concerned but the application of this rule resiilts in a change 
of plan form. In the fictitious incompressible plane, the span 
decreases while the axial distances remain constant, causing the angle 
of sweep to increase. The physical sweep angle A of a given wing at 
a given value of M^ in the compressible plane then corresponds to a 
value of sweep for M^^ = 0, namely, Aq, which is required by the 
present method. Figure 11 shows the variation of with M-^ for 

different values of A. Interpolation in figures 3 and 6 will give the 
corrv^ct velocity ratio at any Mach number for a given wing. 

There are some limitations of this blocking-correction procedure 
which must also be considered. It will be noted that in each case the 
axial distance between the sources and the origin or between the sources 
and sinks approaches zero. Therefore, it has tacitly been assumed that 
the body length or chord must be small compared with the dimensions of 
the tunnel. Should such not be the case,, a small correction must be 
applied to the velocity ratio. The magnitude and nature of this cor- 
rection for a circular wind tunnel are given in reference 1, For the 
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small models normally encountered in high-subsonic-speed testing, this 
correction may usually be neglected. The derivations have also been 
carried out for constant-chord wings, while most models have tapered 
wings. The effect of taper may be taken into account by assuming a 
spanwise variation of g, but this would involve much additional alge- 
braic complexity, particularly for swept wings, which would multiply 
the amount of work in the series summations. Furthermore, only the 
velocity ratios at the centers of the wings have been calculated, while 
the mean along the span should be used for any given wing. For straight 
wings, the additional calculations could easily be carri’ed out using 
equation ( 19 ). For swept wings, however, each position along the span 
is in a different axial plane and would require an additional set of 
equations of greater complexity. The errors involved in neglecting the 
effects of taper and of the spanwise variation of velocity ratio tend 
to compensate each other. 

In conclusion, it should be mentioned that these results may 
easily be applied to another type of testing frequently employed at the 
Cooperative Wind Tunnel. This is the semispan-reflection-plane method 
wherein the model is mounted on the tunnel floor directly without any 
strut support system and vvith the wing extending vertically from the 
floor. It is to be noted that the procedure simply doubles €he effec- 
tive test-section area but does not alter the dimensional ratio 
H/W = l/'| 2 . Therefore, the preceding results are directly applicable 
to a reoriented test section in which A is the center line of the tiuinel 
floor on which the model is mounted, B and E are in the floor corners, 
and C is then in the center line of the tunnel ceiling. 


Southern California Cooperative Wind Tunnel 
California Institute of Technology 
Pasadena, Calif., May 21 , I9U8. 



NACA TN 2372 


23 


APPENDIX A 

CALCULATION OF REMIITOERS 


As mentioned previously, the convergence of the series derived in 
the text of this paper is, in general, quite slow. Thus for practical 
purposes it is necessary to carry the numerical calc-ulations for the 
contributions due to individual images to some reasonable values of n 
and m, say, n = m = n]_. Furthermore, if it can be assumed that n]_ 

is large, certain simplifying assumptions can be applied to the series 
equations, enabling a direct integration in both the m- and n-directions 

1 X 

to be performed from ^ to infinity. The value n^ + inci- 

dentally, comes about from the fact that the individual series summation 
for n^ extends to the outer boundary of that particular image and 
hence, since the point n = m = n]_ is in the center of the image, the 
outer boundary, at which the integration must begin, is at n^_ + 

Two cases are considered, namely, that of the body of revolution 
in the center of the tunnel, and the straight, untapered, finite-span 
wing. It is readily apparent that the effect of sweep rapidly becomes 
negligible as the distance from the primary image becomes large compared 
with the dimensions of a swept wing. Since n^_ is assumed large, the 

remainder terms for swept wings may be considered to be identical to 
those for straight wings. The very rapid convergence of the swept to 
the straight case is immediately apparent in the calculations for indi- 
■vidual images for which m and n are less than n^. A value of seven 
was used for n]_ in all niomerical calculations. 

In all cases the remainders have been broken into separate parts 
in a similar manner. Since nQ_ is assumed to be large, it follows 

that the remainders are independent of the particular location in the 
primaiy image (i.e.. A, B, C, etc.) under consideration. Hence, it is 
sufficient to calculate the results for only one quadrant and simply 
mtiltiply by four. In the actual calculations, an integration was 

carried out from ^ to infinity in both the n- and m-directions 

and the result miltiplied by four. Then the image strips for m = 0 
were integrated from n]_ + i to infinity in the n-direction and multi- 
plied by two. This was repeated for the individual strips for m = 1 
through n^ in the n-direction and the results for each were multiplied 
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by four. The same procedure was carried out for the individual image 
strips in the m-direction. The above may be readily seen in the 
following analyses. 


Body of Revolution 

The velocities at A, B, and C due to the solid blocking effect of 
a body of revolution in the center of a rectangular wind tunnel are 
given by equations ( 8 ), ( 9 ), and (lO), respectively. Using the pro- 
cedure outlined in the preceding paragraphs ^ the results are 
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The remainder terms are the double-integral term and the two braced 


quantities^ 



under the assumptions used, 
results in 


, and are identical for A^ and C 


Performance of the indicated integration 
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Therefore^ the remainder is a function only of rij^ and the dimensional 

ratio of the tunnel H/W and may readily be calculated for any par- 
ticular case« 

Untapered^ Finite-Span Wing 

As mentioned previously^ the effect of sweep upon the remainder 
as n^ becomes large is negligible. Therefore^ the simpler equations 

for Aq = 0 may be integrated and applied to the swept-wing case as 

well. The velocities at A^ and C due to the solid blocking effect 
of a straight^ finite-span wing are given by equations (16)5 (17) ^ and 
(18)^ respectively^ Again applying the assumption that nq is large^ 
these equations may be transformed into 
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The remainder terms again are the double-integral term and the two 
braced quantities, < ^ and \ , and are identical for A, B, 

LJRn LJRni 

and C under the assumptions used. Performance of the indicated inte- 
gration leads to the result 



Since the remainder for the wing case is a function only of n^, H/W, 
and s/W, it may be calculated for any particular case. 
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Figure Two-dimensional wing. 
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Figure Variation of u. /u't, with span. 


H/W = l/\f2 
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Figure ?. 



- Nondimensional velocity- at A due to a line doublet as 
a function of span* 









Angle of sweep, Aq, deg 

Figure 10. -• Variation of ratio of u^ with sweep to without sweep 

with angle of sweep. 






Mach number, M^ 


Figure 11.- Variation of Aq/A with Mach number 
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